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Abstract 

A new theoretical framework is proposed for the description of baryon hght-cone 
distribution amphtudes, based on the observation that their scale dependence to 
leading logarithmic accuracy is described by a completely integrable model. The 
physical interpretation is that one is able to find a new quantum number that 
distinguishes partonic components in the nucleon with different scale dependence. 



1 Introduction 



The notion of baryon distribution amplitudes refers to the valence component 
of the Bethe-Salpeter wave function at small transverse separations and is 
central for the theory of hard exclusive reactions involving baryons [1]. As usual 
for a field theory, extraction of the asymptotic behavior (here: zero transverse 
separation) introduces divergences that can be studied by the renormalization- 
group (RG) method. The distribution amplitude thus becomes a function of 
the three quark momentum fractions Xj and the scale that serves as a UV cutoff 
in the allowed transverse momenta. Solving the corresponding RG equations 
one is led to the expansion 



0(a;„ = 120x1X2X3 E MQl)PN,,ixi) (1) 

N,q \"sWoJ/ 

where the summation goes over all multiplicatively renormalizable operators 
built of three quarks and N derivatives. The polynomials PN,q{xi) and anoma- 
lous dimensions 7^^^^ are obtained by the diagonalization of the mixing matrix 
for the three-quark operators 
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Fig. 1. The spectrum of anomalous dimensions = {l-\-l/Nc)Ei\f + 3/2Cf for the 
baryon distribution ampUtudes with heUcity A = 1/2. The Hues of the largest and 
the smallest eigenvalues for A = 3/2 are indicated by dots for comparison. 



and 4'N,q{Q'i) are the corresponding (nonperturbative) matrix elements. 

As well known, conformal symmetry allows one to resolve the mixing with 
operators containing total derivatives [2]-[9] and requires that the 'eigenfunc- 
tions' Pns corresponding to different values of are mutually orthogonal 
with the weight function 120x1X2X3 that plays the role of the asymptotic wave 
function. In difference to mesons, the conformal symmetry is not sufficient, 
however, to solve the RG equations: for each there exist A" + 1 indepen- 
dent operators with the same conformal spin which mix with each other. This 
mixing produces a nontrivial spectrum of anomalous dimensions, see Fig. 1, 
about which little was known until recently. In what follows, I describe a new 
approach to the scale dependence of baryon distribution amplitude developed 
in [11,12]. The main result is that one is able to identify the second summation 
index q in (1) with an eigenvalue of a certain conserved charge. The physical 
interpretation is that one is able to find a new 'hidden' quantum number that 
distinguishes between partonic components in the proton with different scale 
dependence. 
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2 Hamiltonian formulation of the evolution equations 



In the usual formulation of the evolution equations the conformal invariance 
is not explicit. It can be made manifest, however, if the evolution kernels 
are rewritten in terms of the generators of the SL(2) collinear subgroup of 
conformal transformations. The scale dependence is different for helicity A = 
3/2 and A = 1/2 operators and the corresponding evolution kernels can be 
written in the following compact form [11,12]: 



H- 



3/2 



i<k 



H 



1/2 



i^3/2 - 2 ( 1 + — 



J,fe)-V^(2)J + -Cp, 



+ 



J\2{Jl2 — 1) J23{J23 — 1) 



(2) 
(3) 



Here ip{x) is the logarithmic derivative of the F-function and Jik,i,k = 1,2,3 
are defined in terms of the two-particle Casimir operators of the SL{2, R) 
group 



MJik-l) = Lif, = {Li + Lk)\ (4) 

— * 

with Li being the group generators acting on the i-th quark. In the momentum 
fraction representation (1) the generators take the form [12] 



LkflP{xi) = {xkdk + l)P{xi) , 
Lk,+P{xi) = -XkP{xi) , 

Lk,-P{xi) = {xkdl + 2dk)P{xi). (5) 

Solution of the evolution equations corresponds in this language to solution of 
the Schrodinger equation 



HPn^^{xi)^^n,^P^AXi) (6) 

with 7Ar^q being the anomalous dimensions. The SL{2,R) invariance of the 
evolution equations implies that the generators of conformal transformations 
commute with the 'Hamiltonians' 



[H,L']^[H,L^]^0, (7) 

where = {Li + L2 + L^)^ and = Li a+L^^a + L^ a, so that the polynomials 
PN,q{xi) corresponding to multiplicatively renormalizable operators can be 
chosen simultaneously to be cigcufunctions of and Lq: 
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L^PN,g^{N + S){N + 2)PN,g, 
LoP^,g = {N + 3)Pjv,, , 



(8) 



The third condition in (8) ensures that the operators do not contain overall 
total derivatives. 

Main finding of [11] is that the Hamiltonian H2/2 possesses an additional 
integral of motion (conserved charge): 

Q = '^[Ll2,Ll,] = i{d,-d2){d2-ds){ds-di)xiX2Xs, [H3/2,Q] = 0. (9) 

The evolution equation for baryon distribution functions with maximum heli- 
city is, therefore, completely integrable. The premium is that instead of solving 
a Schrodinger equation with a complicated nonlocal Hamiltonian, it is suffi- 
cient to solve a much simpler equation 



QPN,g{xi) = qPN,,{xi) . (10) 

Once the eigenfunctions are found, the eigenvalues of the Hamiltonian (anoma- 
lous dimensions) are obtained as algebraic functions of N, q. 

It is necessary to add that the Hamiltonian in (2) is known as the Hamiltonian 
describing the so-called XXXs=-i Heisenberg spin magnet. The same Hamil- 
tonian was also encountered in the theory of interacting reggeons in QCD 
[13,14]. 



3 Summary of results: Helicity A = 3/2 distributions 



The equation in (10) cannot be solved exactly, but a wealth of analytic results 
can be obtained by means of the 1/N expansion [15,12]. The general structure 
of the spectrum is illustrated in Fig. 2. It is easy to see that if q is an eigenvalue 
of Q, then —q is also an eigenvalue, whereas the Hamiltonian only depends 
on the absolute value \q\. It follows that all anomalous dimensions arc double 
degenerate except for the lowest ones for each even N, corresponding to the 
solution with q = 0. The corresponding eigenfunctions have a very simple form 

XiX2X3PN,q=o{xi) =xi{l - Xi)C^Jl^{l - 2a;i) + (1 - X2)C^+i(l - 2x2) 
+ xs{l-xs)C^Jl^{l-2xs) (11) 

and the anomalous dimension is equal to 
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Fig. 2. The spectrum of eigenvalues for the conserved charge Q (a) and for the 
hehcity-3/2 Hamiltonian i?3/2 (b), see text. Notation: h = N + 3. 

7iv,,=o = (1 + 1/iVe) [4V^(iV + 3) + - g] + 3/2Cf . (12) 

Furthermore, the eigenvalues of Q he on trajectories (see Fig. 2) corresponding 
to the semi-classicaUy quantized sohton waves [16]. The corresponding trajec- 
tories for the anomalous dimensions have a rather peculiar form. Each of them 
can be considered as a separate partonic component in the nucleon wave func- 
tion, in the same spirit as the leading- twist parton distribution functions in 
deep-inelastic scattering arises from the analytic continuation of the anoma- 
lous dimensions, giving rise to the Altarelli-Parisi splitting function. The 
asymptotic expansions for the charge q and the anomalous dimensions at large 
are available to the order 1/A^^ [15,12] and give very accurate results. The al- 
gebraic structure of the spectrum is very complicated. As an example, I present 
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Fig. 3. The dependence of the energy E on the charge g for iV = 30. The sohd 
curve is calculated using Eq. (13) and the exact values of energy for quantized q are 
shown by crosses. 

the expression for the anomalous dimension jN,q = (1 + ^/Nc)EN,q + 3/2Cf 
as a function of the charge q: [15] 



^iv,g = 2 In 2 - 6 + 67E + 2Re ^(1 + t7j%) + 0(r/-^) , (13) 

k=l 



where rj = (N + 3){N + 2) and S^, k = 1, 2, 3 are the three roots of the cubic 
equation 



2Sl -dk- q/rf. 
Accuracy of (13) is excellent, as illustrated in Fig. 3. 



(14) 



4 Helicity A = 1/2 distributions 



The additional term in spoils integrability but can be considered as a 
small correction for the most part of the spectrum [12]. Its effect on the two 
lowest levels is drastic, however. To illustrate this, consider the flow of energy 
levels for the Hamiltonian H{e) = J2i<k 4^{.Jik) — ip{2) — e I/LI2 + l/-^23 
(cf. (2), (3)) as a function of an auxiliary parameter e, see Fig. 4. It is seen 
that the two lowest levels decouple from the rest of the spectrum and are 
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Fig. 4. The flow of energy eigenvalues for the Hamiltonian H{e) for N = 30. The solid 
and the dash-dotted curves show the parity-even and parity-odd levels, respectively. 
The two vertical dashed lines indicate i?3/2 = H{e = 0) and H1/2 = H{e = 1), 
respectively (up to the color factors). The horizontal dotted line shows position of 
the 'ground state' given by Eq. (12). 

separated from it by a finite mass gap. As shown in [12], this phenomenon can 
be interpreted as binding of the two quarks with opposite helicity and forming 
a scalar "diquark''. The effective Hamiltonian for the low-lying levels can be 
constructed and turns out to be a generalization of the famous Kroning- Penney 
problem for a particle in a 5-function type periodic potential. The value of the 
mass gap between the lowest and the next-to-lowest anomalous dimensions at 
N ^ 00 can be calculated combining the small-e and the large-e expansions 
and is equal to 

A7 = 0.32-(l + l/iVc) (15) 
in agreement with the direct numerical calculations. 



5 Further developments 



The approach developed in [11,12] turns out to be general and is applicable 
to the analysis of all three-parton systems in QCD, albeit with some modifi- 
cations. The extension to three three-gluon operators is straightforward but 
tedious, since elementary fields have a higher conformal spin [17]. Probably 
more interesting from both phenomenological and mathematical point of view 
are the applications to quark-antiquark-gluon systems. The corresponding evo- 
lution equations are integrable in the limit of large number of colors [11] and 
reduce to a different type of integrable models - the so-called open chains. 
A rather detailed analysis of integrable quark-antiquark-gluon operators has 
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been given in [11,18,19] and, most recently, a concrete phenomenological ap- 
plication of this method to the structure function g2{x, Q^) was considered in 
work [20]. 

To summarize, a powerful mathematical framework has been developed for 
the description of the evolution of baryon distribution amplitudes. The math- 
ematical structure of evolution equations is very elegant and reveals certain 
qualitative features of the distribution amplitudes. A lot of analytical results 
is obtained, in different limits. The formalism is general and was applied al- 
ready to the other existing three-parton distributions. In a more general con- 
text, the integrability of evolution equations reveals an additional symmetry of 
QCD and its close relation to exactly solvable statistical models. Remarkably 
enough, the same symmetry has been observed in the studies of the Regge 
asymptotics of three-gluon distributions. All these features are not seen at the 
level of QCD Lagrangian and their origin has to be understood better. 
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